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1 Introduction 

Statistical mechanical properties of random surfaces as well as randomly 
branched polymers can be analyzed within the framework of large N expan- 
sion. In the same manner in which matrix models in their double scaling limit 
[1,2,3,4] provide representations of dynamically triangulated random surfaces 
summed on different topologies, O(N) symmetric vector models represent 
discretized branched polymers in this limit, where N — > oo and the coupling 
constant g — > g c in a correlated manner [5,6]. The surfaces in the case of ma- 
trix models, and the randomly branched polymers in the case of vector models 
are classified by the different topologies of their Feynman graphs and thus by 
powers of 1/N. Though matrix theories attract most attention, a detailed un- 
derstanding of these theories exists only for dimensions d < 1. On the other 
hand, in many cases, the O(N) vector models can be successfully studied also 
in dimensions d > 1 [6], and thus, provide us with intuition for the search for 
a possible description of quantum field theory in terms of extended objects 
in four dimensions, which is a long lasting problem in elementary particle 
theory. 

The double scaling limit in O(N) vector quantum field theories reveals an 
interesting phase structure beyond N — > oo limit. In particular, though the 
N — > oo multicritical structure of these models is generally well understood, 
there are certain cases where it is still unclear which of the features survives 
at finite N, and to what extent. One such problem is the multicritical be- 
havior of O(N) models at critical dimensions [7]. Here, one finds that in the 
N — > oo limit, there exists a non-trivial UV fixed point, scale invariance is 
spontaneously broken, and the one parameter family of ground states con- 
tains a massive vector and a massless bound state, a Goldstone boson-dilaton. 
However, since it is unclear whether this structure is likely to survive for fi- 
nite N [8], one would like to know whether it is possible to construct a local 
field theory of a massless dilaton via the double scaling limit, where all orders 
in 1/N contribute. The double scaling limit is viewed as the limit at which 
the attraction between the O(N) vector quanta reaches a value at g — > g c , at 
which a massless bound state is formed in the N — > oo limit, while the mass of 
the vector particle stays finite. In this limit, powers of 1/N are compensated 
by IR singularities and thus all orders in 1/N contribute. 

In section 2 the double scaling limit for simple integrals and quantum me- 
chanics is recalled, introducing a formalism which will be useful for field 
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theory examples. 

In section 3 the special case of field theory in dimension two is discussed, 
slightly generalizing previous established results. 

In higher dimensions a new phenomenon arises: the possibility of a spon- 
taneous breaking of the O(N) symmetry of the model, associated to the 
Goldstone phenomenon. 

Before discussing a possible double scaling limit, the critical and multi- 
critical points of the O(N) vector model are reexamined in section 4. In 
particular, a certain sign ambiguity that appears in the expansion of the gap 
equation is noted, and related to the existence of the IR fixed point in dimen- 
sions 2 < d < 4. In section 5 we discuss the subtleties and conditions for the 
existence of an O(N) singlet massless bound state along with a small mass 
O(N) vector particle excitation. It is pointed out that the correct massless 
effective field theory is obtained after the massive O(N) scalar is integrated 
out. Section 6 is devoted to the double scaling limit with a particular em- 
phasis on this limit in theories at their critical dimensions. In section 7 the 
main conclusions are summarized. 

2 Double scaling limit: simple integrals and quantum mechanics 

The double scaling limit [1,2,3,4] of the vector model has already been investi- 
gated [6], for dimensions d < 1 and for the simple (4> 2 ) 2 field theory. We first 
recall results obtained in d = and d — 1 dimensions, dimensions in which 
the matrix models have equally been solved. We however introduce a general 
method, not required here, but useful in the general field theory examples. 

2.1 The zero dimensional example 

Let us first recall the zero dimensional example. The partition function Z is 
given by 



The simplest method for discussing the large iV limit is of course to integrate 
over angular variables (see appendix Al). Instead we introduce two new 
variables A, p and use the identity 



The integral over A is really a Fourier representation of a 5-function and thus 
the contour of integration runs parallel to the imaginary axis. The identity 
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(2.1) transforms the action into a quadratic form in <p. Hence the integration 
over <p can be performed and the dependence in N becomes explicit 

e z oc J dpd\exp{-N [-\\p + V{p) + 5 In A]} . 

The large N limit is obtained by steepest descent. The saddle point is given 
by 

V'(p) = l\, P = i/x. 

The leading contribution to Z is proportional to N and obtained by replacing 
A, p by the saddle point value. The leading correction is obtained by expand- 
ing A, p around the saddle point and performing the gaussian integration. It 
involves the determinant D of the matrix M of second derivatives 

M=("H" 2 y^ p) ). D = detM = -i(V"(p)/A 2 + ±). 

In the generic situation the resulting contribution to Z is — -| lnD. However 
if the determinant D vanishes the leading order integral is no longer gaussian, 
at least for the degree of freedom which corresponds to the eigenvector with 
vanishing eigenvalue. The condition of vanishing of the determinant also 
implies that two solutions of the saddle point equation coincide and thus 
corresponds to a surface in the space of the coefficients of the potential V 
where the partition function is singular (see appendix Al for details). 

To examine the corrections to the leading large N behaviour it remains 
however possible to integrate over one of the variables by steepest descent. 
At leading order this corresponds to solving the saddle point equation for one 
of the variables, the other being fixed. Here it is convenient to eliminate A 
by the equation A = 1/p. One finds 

e z oc J dp exp [-N(V(p) - \ lnp) + 0(1)] . 

In the leading term we obviously recover the result of the angular integration 
with p = <f 2 . For N large the leading contribution arises from the leading 
term in the expansion of W(p) = V(p) — \ lnp near the saddle point: 

w{p)-w{ Ps )^x { w^{p s ){p-p s Y. 

The integer n characterizes the nature of the critical point. Adding relevant 
perturbations 5kV of parameters Vk to the critical potential 

SkV = v k (p - p s ) k , l<k<n-2 
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(the term k = n — 1 can always be eliminated by a shift of p) we find the 
partition function at leading order for N large in the scaling region: 



n-2 



,Z({u k }) 



U k Z 



k=l 



oc J dz exp ^— z n — 
where z oc A^ 1 / n (p — p s ) and 

u k oc N 1 ~ k / n Vf t 

is held fixed. 



Quantum mechanics 

The method we have used above immediately generalizes to quantum me- 
chanics, although a simpler method involves the Schrodinger equation. We 
consider the euclidean action 



S(f) = N / dt \{d t ${t)) 2 + V(0 2 ) 



(2.2) 



Note the unusual field normalization, the factor iV in front of the action 
simplifying all expressions in the large iV limit. 

To explore the large iV limit one has to take the scalar function </> 2 , which 
self-averages, as a dynamical variable [9] . At each time t we thus perform the 
transformation (2.1). One introduces two paths p(t), X(t) and writes 



exp 



N J dtV{(j) 2 ) 
c f [dp(t) d\{t)] 



exp <-N dt 



) + V(p) 



(2.3) 



The integral over the path <f>(t) is then gaussian and can be performed. One 
finds 

z 



with 



J dp(t)dX(t) exp [-^(A, p)] 
S eS = N J dt [-|Ap + y(p) + |trln(-d f 2 + A(t))] . 



Again, in the large limit the path integral can be calculated by steepest 
descent. The saddle points are constant paths solution of 



V'{p) = \X, 



1 



duo 



1 



2tt J uj 2 + X 2VX ' 
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where ui is the Fourier energy variable conjugate to t. Again a critical point 
is defined by the property that at least two solutions to the saddle point 
equations coalesce. This happens when the determinant of the matrix of 
second derivatives of the equations vanishes: 



\ 2 4tt / (o> 2 +A) 2 / 

The leading correction to the saddle point contribution is given by a gaussian 
integration. The result involves determinant of the operator second derivative 
of S e ff. By Fourier transforming time the operator becomes a tensor product 
of 2 x 2 matrices with determinant D(u>) 



D{u) = det ( ' T ! f 2 d , 



V"(p) 

1 jl r 

2 4tt J (lj' 2 +\)[(lo-uj') 2 +X] 

Thus, the criticality condition is equivalent to D(0) = 0. When the criticality 
condition is satisfied, the leading correction is no longer given by steepest 
descent. Again, since at most one mode can be critical, we can integrate 
over one of the path by steepest descent, which means solving the saddle 
point equation for one function, the other being fixed. While the p-equation 
remains local, the A is now non-local, involving the diagonal matrix element 
of the inverse of the differential operator — d 2 + X(t). We shall see in next 
section how this problem can be overcome in general. A special feature of 
quantum mechanics, however, is that the determinant can be calculated, after 
a simple change of variable. We set 

X(t) = s(t) + s 2 (t), 

in such a way that the second order differential operator factorizes 

-d 2 + X(t) = -{d t + s(t)) (d t - s(t)). (2.4) 

The determinant of a first order differential operator can be calculated 

lndet(-d t 2 + X(t)) = J dts(t). 

The jacobian of the transformation (2.4) contributes at higher order in 1/N 
and can be neglected. Therefore the effective action becomes 

S eS = N J dt + s 2 )p + V(p) + ±s(t)] 

= N Jdt [-\ps 2 + l s (p + l) + V(p)]. 
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We can now replace s by the solution of a local saddle point equation: 



-sp+±(p+l) = 0, 

and find 



S eS = N dt 



We recognize the action for the large N potential at zero angular momentum 
in the radial coordinate p(t) = (p 2 (t). Critical points then are characterized 
by the behaviour of the potential W(p) 

w{ P ) = v{ P ) + i- , 

near the saddle point p s 

w( P ) - w(p 8 ) ~ w( n \ Ps ) {p ~ p ° )n . 

nl 

At critical points the ground state energy, after subtraction of the classical 
term which is linear in N, has a non-analytic contribution. To eliminate N 
from the action we set 

t tN (n-2)/(n+2)^ p{f) _ ^ ^ N ~V^+2) Z ^ 

We conclude that the leading correction to the energy levels is proportional 
to N~( n ~ 2 ^( n+2 \ Note also that the scaling of time implies that higher order 
time derivatives would be irrelevant, an observation which can be used more 
directly to expand the determinant in local terms, and will be important in 
next section. 

If we add relevant corrections to the potential 

SkV = v k (p-p s ) k , l<k<n-2, 

the coefficients Vk must scale like 
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3 The 2D V((p 2 ) field theory in the double scaling limit 



In the first part we recall the results concerning the O(N) symmetric V{4> 2 ) 

— * 

field theory, where <p is iV-component field, in the large N limit in dimension 
two because phase transitions occur in higher dimensions, a problem which 
has to be considered separately. The action is: 

S($) = N J d 2 *{± [«9 M #r)] 2 + y (3.1) 

where an implicit cut-off A is always assumed below. Whenever the explicit 
dependence in the cut-off will be relevant we shall assume a Pauli-Villars's 
type regularization, i.e. the replacement in action (3.1) of —(f)d 2 (j) by 

-$d 2 D(-d 2 /A 2 )$, (3.2) 

where D(z) is a positive non-vanishing polynomial with D(0) = 1. 

As before one introduces two fields p(x) and X(x) and uses the identity 
(2.3). The effective action is then: 

S eS = N J d 2 x [V(p) - \\p] + §ATtrm(-A + A). (3.3) 

Again for N large we evaluate the integral by steepest descent. Since the 
saddle point value A is the ^"-field mass squared, we set in general A = m 2 . 
With this notation the two equations for the saddle point m 2 , p s = ((p 2 ) are: 



V'{ Ps ) = \m 2 , (3.4a) 



(27r) 2 J k 2 + m 2 
where we have used a short-cut notation 

1 f A d 2 k = J_ f d 2 k 2 

(2tt) 2 J k 2 + m 2 ~ (2tv) 2 J D(k 2 /A 2 )k 2 + m 2 ~ l[m >' 

For m<A one finds 

Si(m 2 ) = — ln(A/m) + — ln(87rC) + 0(m 2 /A 2 ) 



»=r ds fe" e(i 



As we have discussed in the case of quantum mechanics a critical point is 
characterized by the vanishing at zero momentum of the determinant of sec- 
ond derivatives of the action at the saddle point. The mass-matrix has then 
a zero eigenvalue which, in field theory, corresponds to the appearance of a 

— * 

new massless excitation other than <p. In order to obtain the effective action 
for this scalar massless mode we must integrate over one of the fields [10]. In 
the field theory case the resulting effective action can no longer be written 
in local form. To discuss the order of the critical point, however, we only 
need the action for space independent fields, and thus for example we can 
eliminate A using the A saddle point equation. 
The effective p potential W(p) then reads 

/-Hp) f) 

W{p) = V{p)-\J dX' X'—Bi(X'), (3.5) 

where at leading order for A large 

A(p) = 8nCA 2 e~ Anp . 

The second term in Eq. (3.5) in fact is the kinetic energy contribution to 
the ground state free energy as can be viewed in an Hartree-Fock variational 
calculation that becomes exact in the limit of N — > oo (see for example Ref. 
[11] ) . The expression for the effective action in Eq. (3.5) is correct for any 
d and will be used also in section 6. 
Here we have: 



W(p) = V(p) + CA 2 e" 4 ^ = V(p) + ±m 2 q -*<p-p°) . 
A multicritical point is defined by the condition 



W(p)-W(p s ) = 0((p-p s ) n ) (3.6). 

This yields the conditions: 

V ik) (p s ) = i(-4 7 r) fc - 1 m 2 for 1 < k < n- 1. 

Note that the coefficients V^ k \p s ) are the coupling constants renormalized at 
leading order for N large. If V(p) is a polynomial of degree n — 1 (the minimal 
polynomial model) the multicritical condition in Eq. (3.6) and Eq. (3.46) 
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determines the critical values of n — 1 renormalized coupling constants as well 
as p s in terms of m 2 . 

When the fields are space-dependent it is simpler to eliminate p instead, 
because the corresponding field equation: 

V(p{x)) = \X{x). (3.7) 

is local. This equation can be solved by expanding p(x) — p s in a power series 
in X(x) — m 2 : 

p(x) - Ps = ^L-y (X(x) - m 2 ) + O ((A - m 2 ) 2 ) . (3.8) 

The resulting action for the field X(x) remains non-local but because, as we 
shall see, adding powers of A as well as adding derivatives make terms less 
relevant, only the few first terms of a local expansion of the effective action 
will be important. 

If in the local expansion of the determinant we keep only the two first terms 
we obtain an action containing at leading order a kinetic term proportional 
to (<9 M A) 2 and the interaction (X(x) — m 2 ) n : 

S eff (A) ~ N J d 2 x ^-L^A) 2 + ±S n (X(x) -m 2 y 

where the neglected terms are of order (A — m 2 ) n+1 , A<9 4 A, and X 2 d 2 X and 

S n = W^(p s )[2V"(p 8 )]- n = W^(p s )(-4nm 2 )- n . 

Moreover we note that together with the cut-off A, m now also acts as a 
cut-off in the local expansion. 

To eliminate the N dependence in the action we have, as in the example of 
quantum mechanics, to rescale both the field A — m 2 and coordinates: 

X(x) - m 2 = v / 487m 2 A^- 1/2 v?(x) , x i-> N (n ~ 2) / 4 x . (3.9) 

We find 



SeffM- J d 2 x [l(d„<p) 2 + ±g n <p n ] 



In the minimal model, where the polynomial V(p) has exactly degree n — 1, 
we find g n = 6(487r)( n - 2 )/ 2 m 2 . 
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As anticipated we observe that derivatives and powers of tp are affected by 
negative powers of N, justifying a local expansion. However we also note that 
the cut-offs (A or the mass m) are now also multiplied by jV~( n_2 )/ 4 . Therefore 
the large N limit also becomes a large cut-off limit. 

Double scaling limit. The existence of a double scaling limit relies on the 
existence of IR singularities due to the massless or small mass bound state 
which can compensate the 1/N factors appearing in the large N perturbation 
theory. 

We now add to the action relevant perturbations: 

S k V = v k (p(x) - p s ) k , \<k<n-2. 
Namely, adding to the A action a sum of terms proportional to J d 2 x(X—m 2 ) k : 

2\k 



S k S eS (X)=NS k J d 2 x(X-m 2 ) 



where the coefficients S k are functions of the coefficients v k . After the rescal- 
ing (3.9) these terms become 

4^eff(v) = H^ iv(n " fc)/2 / dW(*0 l<k<n-2 

However, unlike quantum mechanics, it is not sufficient to scale the coeffi- 
cients g k with the power N^ h ~ n ^ 2 to obtain a finite scaling limit. Indeed 
perturbation theory is affected by UV divergences, and we have just noticed 
that the cut-off diverges with N. In two dimensions the nature of divergences 
is very simple: it is entirely due to the self-contractions of the interaction 
terms and only one divergent integral appears: 

where p is the small mass of the bound state, required as an IR cut-off to 
define perturbatively the double scaling limit. We can then extract the N 
dependence 

(<p 2 (x)) = — (n - 2) In N + 0(1). 

Therefore the coefficients S k have also to cancel these UV divergences, and 
therefore have a logarithmic dependence in N superposed to the natural power 
obtained from power counting arguments. In general for any potential U (<p) 

2k' 



U(<p) =: U(<p) : + 



^ 2 k k\ ^ ^ (<9v?) 



: U(<p) : , 
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where : U(ip) : is the potential from which self-contractions have been sub- 
tracted (it has been normal-ordered). For example for n = 3 

ip 3 (x) =: ip 3 (x) : +3(ip 2 )ip(x), 

and thus the double scaling limit is obtained with the behaviour 

Ngi H — — IniV and N 1/2 g 2 fixed . 

167T 

For another example n = 4 



#iiV 3/2 and Ng 2 + —InN fixed 

87T 



4 The V{4> 2 ) in the large N limit: phase transitions 



In higher dimensions something new happens: the possibility of phase transi- 
tions associated with spontaneous breaking the O(N) symmetry. In the first 
part we thus study the O(N) symmetric V((p 2 ) field theory, in the large iV 
limit to explore the possible phase transitions and identify the corresponding 
multicritical points. The action is: 

S($) = N J d d x | \ [dj(x)] 2 + V (V 2 ) } , (4.1) 



where, as above (Eqs.(3.1)-(3.2)), an implicit cut-off A is always assumed 
below. 

The identity (2.3) transforms the action into a quadratic form in <p and 
therefore the integration over <p can be performed. It is convenient however 
here to integrate only over N — 1 components, to keep a component of the 
vector field, which we denote cr, in the action. The effective action is then: 



S eS = N d d x 



\ (<V) 2 + V(p) + \\ (a 2 - p)] + \(N - l)trln(-A + A). 

(4.2) 
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4-1 The saddle point equations: the O(N) critical point 

Let us then write the saddle point equations for a general potential V. At 
high temperature a = and A is the 0-field mass squared. We thus set in 
general A = m 2 . With this notation the three saddle point equations are: 



mV = 0, (4.3a) 



V'(p) = \m 2 , (4.36) 



2 1 f A d d k 



with the notation of section 3 



1 f A d d k If d d k _ 2 

(27r) d J k 2 + m 2 ~ (2n) d J D(k 2 /A 2 )k 2 + m 2 ~ ll ™ h 

In the ordered phase a ^ and thus m vanishes. Equation (4.3c) has a 
solution only for p > p c , 

" A d d k 



i r d a k . 



(2n) 

Equation (4.36) which reduces to V'(p) = then yields the critical tempera- 
ture. Setting V(p) = U(p) + \rp, we find 

r c = -2U'(p c ). 

To find the magnetization critical exponent (3 we need the relation between 
the r and p near the critical point. 

In the disordered phase, a = 0, equation (4.3c) relates p to the 0-field mass 
m. It can be rewritten 

p = B 1 (m 2 ) = A d " 2 F(m 2 /A 2 ), (4.4) 

where 

, . 2 f k dk , . 

^ = (An) d / 2 T{d/2) J k 2 D(k 2 )+z ' ( 

The function -F(^) can be written in an asymptotic expansion 

oo oo 

= z d/2 - 1 ^b n z n + ^c n z n . (4.6) 

n=0 n=0 
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The non-analytic part can be extracted from the representation 

d/2-1 poo poo 

F(z) = , ,, w w/0 / dx dyy- d / 2 x d / 2 - 1 exp\--D(zx/y)-?-\, 
KI 2r(rf/2)(4yr) d / 2 J J yy F l 2 v 1 yj 2J' 

(4.7) 

which gives e.g. in the case of D(k 2 ) = 1 + <ii/c 2 the asymptotic expansion for 
the non-analytic part: 

F (z) - r(1 - d/2) (SF5^ o „i r(f +n) ' (48) 

and the analytic part is obtained from -F(z) — F NA (z). For m <C A, p ap- 
proaches p c , 

p c = S 1 (0) = A d - 2 F(0), 

and the relation becomes: 

p - Pc = -K(d)m d - 2 + a(d)m 2 A d - 4 + O (m d A" 2 ) + O (m 4 A d " 6 ) . (4.9) 

For 2 < d < 4 (the situation we shall assume below except when stated oth- 
erwise) the O (m d A -2 ) from the non-analytic part dominates the corrections 
to the leading part of this expression. For d = 4 instead 

p — p c = g^rn 2 (lnm/A + const.) , 

and for d > 4 the analytic contribution dominates and 

p - p c ~ a(d)m 2 A d - 4 . 

The constant if (<i) is universal: 

, 7 * , 1 /" d d A; r(l-d/2) , 

The constant a(d) instead depends on the cut-off procedure, and is given by 



^-^vkpj^k-wm 1 (411) 



Let us also define for later purpose the function 
1 f d d k 



B 2 (p-m 2 ) 



(27r) d J [k 2 D(k 2 /A 2 ) + m 2 } [(p - k) 2 D((p - k) 2 / A 2 ) + m 2 ] ' 

(4.12) 
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which is up to sign the second derivative with respect to A (a?) of the tr In 
term in the effective action. Then for m small 

^™ 2 ) = -H^(™ 2 ) = -5^(<>-<>d 

= \{d- 2)K(d)m d - 4 - a(d)A d - 4 + O (m^A" 2 , m 2 A d " 6 ) . 

(4.13) 

Critical point. In a generic situation V"(p c ) = U"(p c ) does not vanish. We 
thus find in the low temperature phase 

t = r _ rc ^_2U"(p c )(p-p c ) (3=\. (4.14) 

This is the case of an ordinary critical point. Stability implies V"(p c ) > so 
that t < 0. 

At high temperature, in the disordered phase, the </>-field mass m is given 
by 2U'(p) + r = m 2 and thus, using (4.9), at leading order 

t ~ 2U"(p c )K(d)m d - 2 , 

in agreement with the result of the normal critical point. Of course the 
simplest realization of this situation is to take V(p) quadratic, and we recover 
the (4> 2 ) 2 field theory. 

The sign of the constant a(d). A comment concerning the non-universal 
constant a(d) defined in (4.9) is here in order because, while its absolute value 
is irrelevant, its sign plays a role in the discussion of multicritical points. 
Actually this sign is already relevant to the RG properties of the large N 
limit of simple scalar field theories. In a V(4> 2 ) = ^(p 2 + ^{4> 2 ) 2 theory (in 
Eq. (4.1)) it is easy to verify that a(d) is related to the second coefficient 
of the large RG /^-function. If we call A^A = gA A ~ d the bare coupling 
constant, we indeed find [9]: 

(3(g) = -(4 - d)g + 1(4 - d)Na(d)g 2 + O (1/N) . (4.15) 

It is generally assumed that a(d) > 0. Indeed, this is what is found near 
four dimensions in all regularizations. Then there exists an IR fixed point, 
non-trivial zero of the /9-function. For the simplest Pauli-Villars's type regu- 
larization we have D(z) > 1 and thus a(d) is finite and positive in dimensions 
2 < d < 4, but this is not a universal feature. Even in case of simple lattice 
regularizations it has been shown [12] that in d = 3 the sign is arbitrary. We 
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illustrate the ambiguity in the sign of a(d) in 2 < d < 4 by explicit examples 
in appendix A2. However, if a(d) is negative the large N RG has a problem, 
since the coupling flows in the IR limit to large values where the large N 
expansion is no longer reliable. It is not known whether this signals a real 
physical problem, or is just an artifact of the large N limit. 

Another way of stating the problem is to examine directly the relation 
between bare and renormalized coupling constant. Calling g T m 4 ~ d the renor- 
malized 4-point function at zero momentum, we find 

4-d _ iv y i a 

m g r - 1 + A 4-d gNB2 ( . m 2y 6 - ^ ib J 
In the limit m< A the relation can be written 



We see that when a(d) < the renormalized IR fixed point value cannot be 
reached by varying g > for any finite value of m/A. 

4-2 Multicritical points 

A new situation arises if we can adjust a parameter of the potential in such 
a way that U"(p c ) = 0. This can be achieved only if the potential V is at 
least cubic. We then expect a tricritical behavior. Higher critical points can 
be obtained when more derivatives vanish. We shall examine the general case 
though, from the point of view of real phase transitions, higher order critical 
points are not interesting because d > 2 for continuous symmetries and mean- 
field behavior is then obtained for d > 3. The analysis will however be useful 
in the study of double scaling limit. 

Assuming that the first non- vanishing derivative is U^ n \p c ), we expand 
further equation (4.36). In the ordered low temperature phase we now find 

t= -J^^ uin) MiP-Pc) n -\ =► P=^TTy (4-18) 

which leads to the exponent (3 expected in the mean field approximation for 
such a multicritical point. We have in addition the condition U^ n \p c ) > 0. 
In the high temperature phase instead 

m 2 = t + (-I) 71 ' 1 2 _ U^ n \ Pc )K n - l {d)m^- l ^ d -^. (4.19) 
(n — 1). 
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For d > 2n/(n — 1) we find a simple mean field behavior, as expected since 
we are above the upper-critical dimension . 

For d < 2nj (n — 1) we find a peculiar phenomenon, the term in the r.h.s. is 
always dominant, but depending on the parity of n the equation has solutions 
for t > or t < 0. For n even, t is positive and we find 

•»«*". "= (»- !)'(„- 2) - (420) 

which is a non mean-field behavior below the critical dimension. However for 
n odd (this includes the tricritical point) t must be negative, in such a way 
that we have now two competing solutions at low temperature. We have to 
find out which one is stable. We shall verify below that only the ordered phase 
is stable, so that the correlation length of the 0-field in the high temperature 
phase always remains finite. Although these dimensions do not correspond to 
physical situations because d < 3 the result is peculiar and inconsistent with 
the e-expansion. 

For d = 2n/(n — 1) we find a mean field behavior without logarithmic 
corrections, provided one condition is met: 

U^Mir- 1 (2n/(n -!))<!, K(3) = 1/(4tt). (4.21) 



(n- 1)! 



We examine, as an example, in more details the tricritical point below. We 
will see that the special point 

■U^ipjK"- 1 (2n/(n - 1)) = 1 , (4.22) 



(n-1) 

has several peculiarities [7]. In what follows we call fi c this special value of 
U^(p c ). 

Discussion. In the mean field approximation the function U(p) oc p n is not 
bounded from below for n odd, however p = is the minimum because by 
definition p > 0. Here instead we are in the situation where U (p) ~ (p — p c ) n 
but p c is positive. Thus this extremum of the potential is likely to be unstable 
for n odd. To check the global stability requires further work. The question is 
whether such multicritical points can be studied by the large iV limit method. 

Another point to notice concerns renormalization group: The n = 2 ex- 
ample is peculiar in the sense that the large limit exhibits a non-trivial 
IR fixed point. For higher values of n no coupling renormalization arises in 
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the large N limit and the IR fixed point remains pseudo-gaussian. We are 
in a situation quite similar to usual perturbation theory, the (3 function can 
only be calculated perturbatively in 1/N and the IR fixed point is outside 
the perturbative regime. 

4-3 Local stability and the mass matrix 

The matrix of the general second partial derivatives of the effective action is: 

(p 2 + m 2 a \ 

V"(p) -\ , (4.23) 

-\ -lB 2 (p;m 2 )J 

where B2(p;m 2 ) is defined in (4.12). 

We are in position to study the local stability of the critical points. Since 
the integration contour for A = m 2 should be parallel to the imaginary axis, 
a necessary condition for stability is that the determinant remains negative. 

The disordered phase. Then a = and thus we have only to study the 
2x2 matrix M of the p, m 2 subspace. Its determinant must remain negative 
which implies 

detM<0 2V"(p)B 2 (p;m 2 ) + 1 > 0. (4.24) 

For Pauli-Villars's type regularization the function B 2 (p:,m 2 ) is decreasing 
so that this condition is implied by the condition at zero momentum 

detM<0 <= 2V"(p)B 2 (0;m 2 ) + 1 > 0. 

For m small we use Eq. (4.13) and at leading order the condition becomes: 

K{d) (d - 2)m d - 4 V"(p) + 1 > . 

This condition is satisfied by a normal critical point since V"(p c ) > 0. For a 
multicritical point, and taking into account equation (4.9) we find: 

(-l) n ^ ~ 2 2 y K n - l (d)m n{d -V- d V( n) (p c ) + 1 > . (4.25) 

We obtain a result consistent with our previous analysis: For n even it is 
always satisfied. For n odd it is always satisfied above the critical dimension 
and never below. At the upper-critical dimension we find a condition on 
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the value of V^ n \p c ) which we recognize to be identical to condition (4.21) 
because then 2/(n — 1) = d — 2. 

The ordered phase. Now m 2 = and the determinant A of the complete 
matrix is: 

-A>0 2V"(p)B 2 (0;p 2 )p 2 + p 2 +4V"(p)a 2 > 0. (4.26) 

We recognize a sum of positive quantities, and the condition is always sat- 
isfied. Therefore in the case where there is a competition with a disordered 
saddle point only the ordered one can be stable. 

5 The scalar bound state 

In this section we study the limit of stability in the disordered phase (a = 0) . 
This is a problem which only arises when n is odd, the first case being provided 
by the tricritical point. The mass-matrix has then a zero eigenvalue which 
corresponds to the appearance of a new massless excitation other than </>. Let 
us denote by M the p, m 2 2 x 2 submatrix. Then 

detM = 2V"(p)B 2 (0;m 2 ) + 1 = 0. 

In the two-space the corresponding eigenvector has components (|, V"(p)). 

5.1 The small mass m region 

In the small m limit the equation can be rewritten in terms of the constant 
K(d) defined in (4.10): 

K(d)(d-2)m d - 4 V"(p) + 1 = 0. (5.1) 

Equation (5.1) tells us that V"(p) must be small. We are thus close to a 
multicritical point. Using the result of the stability analysis we obtain 

(-l)"" 1 ^ ~ ^ ""^jm^-^VW (p c ) = 1 . (5.2) 

We immediately notice that this equation has solutions only for n(d — 2) = 
d, i.e. at the critical dimension. The compatibility then fixes the value of 
V^(p c ). We again find the point (4.22), V^(p c ) = fi c . If we take into 
account the leading correction to the small m behavior we find instead: 
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This means that when a(d) > there exists a small region V^ n \p c ) > fl c 
where the vector field is massive with a small mass m and the bound-state 
massless. The value fi c is a fixed point value. 

The scalar field at small mass. We want to extend the analysis to a situation 
where the scalar field has a small but non-vanishing mass M and m is still 
small. The goal is in particular to explore the neighborhood of the special 
point (4.22). Then the vanishing of the determinant of M implies 

l + 2V"(p)B 2 (iM;m 2 ) = 0. (5.4) 

Because M and m are small, this equation still implies that p is close to a 
point p c where V"(p) vanishes. Since reality imposes M < 2m, it is easy to 
verify that this equation has also solutions for only the critical dimension. 
Then 

V^(p c )f(m/M) = n c , (5.5) 

where we have set: 

/(*)= / dx[l + (x 2 -l)/(4z 2 )] d/2 -\ \<z. (5.6) 
Jo 

In three dimensions it reduces to: 

'2z + l 



f(z) = zln 



2z 



f(z) is a decreasing function which diverges for z = \ because d < 3. Thus 
we find solutions in the whole region < V ^ (p c ) < fi c , i.e. when the multi- 
critical point is locally stable. 

Let us calculate the propagator near the pole. We find the matrix A 



2 

G2 



N dB 2 (p;m 2 ] 



dp 2 

where we have set 



p2 = _M 2 



1 / 1 G 

p 2 + m 2 \G G 2 ' 



(5.7) 



2( -Kr-^) 

(n-2)! 

For m/M fixed the residue goes to zero with m as m d ~ 2 because the deriva- 
tive of B is of the order of m d ~ 6 . Thus the bound-state decouples on the 
multicritical line. 
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5.2 The scalar massless excitation: general situation 

Up to now we have explored only the case where both the scalar field and the 
vector field propagate. Let us now relax the latter condition, and examine 
what happens when m is no longer small. The condition M = then reads 

2V"{p s )B 2 {0; m 2 ) + 1 = 

together with 

m 2 = 2V'(p s ), p s = B 1 (m 2 ). (5.8) 

An obvious remark is: there exist solutions only for V"(p s ) < 0, and therefore 
the ordinary critical line can never be approached. In terms of the function 
F(z) defined by equation (4.4) the equations can be rewritten 

p s = A d - 2 F(z), z = 2V'( Ps )A- 2 , 2K d - A V"{ Ps )F\z) = 1 . 

The function F(z) in Pauli-Villars's regularization is a decreasing function. 
In the same way —F'{z) is a positive decreasing function. 

The third equation is the condition for the two curves corresponding to the 
two first ones become tangent. For any value of z we can find potentials and 
thus solutions. Let us call z s such a value and specialize to cubic potentials. 
Then 

p s = K d - 2 F(z s ) , 

V(p) = V'(p s )(p - Ps ) + \V"{ Ps ){p - Ps ) 2 + iV^(p s )(p - p s ) 3 , (5.9) 

which yields a two parameter family of solutions. For z small we see that for 
d < 4 the potential becomes proportional to (p — p c ) 3 . 

6 Stability and double scaling limit 

In order to discuss in more details the stability issue and the double scaling 
limit we now construct the effective action for the scalar bound state. We 
consider first only the massless case. We only need the action in the IR limit, 
and in this limit we can integrate out the vector field and the second massive 
eigenmode. 

Integration over the massive modes. As we have already explained in section 
3 we can integrate over one of the fields, the second being fixed, and we need 
only the result at leading order. Therefore we replace in the functional integral 



e z 



= J[dpd\]exp -ftrln(-<9 2 + A) + N J d d x(-V(p) + \p\) 



(6.1) 
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one of the fields by the solution of the field equation. It is useful to discuss 
the effective potential of the massless mode first. This requires calculating 
the action only for constant fields it is then simpler to eliminate A. We 
assume in this section that m is small (the vector propagates). For A < A 
the A-equation reads (d < 4) 

p-p c =-K(d)X^/ 2 . (6.2) 

It follows that the resulting potential W(p), obtained from Eq. (3.5) is 

WM = v(p) + wmw^ (p < ~ p)J/< " 2 ' (6 ' 3) 

In the sense of the double scaling limit the criticality condition is 

w{p) = o({p- Ps T). 

It follows 

y<*>(p a ) = - 1 jri-*(d) Efc ~ d,{d - m d-k(d-2) 1 < k < n _ 1 

yi ! 2 v ; r(-2/(d-2)) - - 

For the potential V of minimal degree we find 

W {p) ~ 9 l T g 1 -"(d) r ^ ~ d /} d " ^ m d -< d - 2 \p - p s ) n . 

\I J ) 2n! V ) r(_ 2 /(d_2)) ^ ^ ^ 



T/ie double scaling limit. We recall here that quite generally one verifies 
that a non-trivial double scaling limit may exist only if the resulting field 
theory of the massless mode is super-renormalizable, i.e. below its upper- 
critical dimension d = 2n/(n — 2), because perturbation theory has to be IR 
divergent. Equivalently, to eliminate iV from the critical theory, one has to 
rescale 

p — p s (x N~ 29 tf , x^xN^ n - 2)9 with 1/9 = 2n-d{n-2), 

where 6 has to be positive. 

We now specialize to dimension three, since d < 3 has already been exam- 
ined, and the expressions above are valid only for d < 4. The normal critical 
point (n = 3), which leads to a </? 3 field theory, and can be obtained for a 
quadratic potential V(p) (the {4> 2 ) 2 ) has been discussed elsewhere [6]. We 
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thus concentrate on the next critical point n = 4 where the minimal potential 
has degree three. 

The d = 3 tricritical point. The potential W(p) then becomes 

W(p) = V(p) + ^(p c -pf. (6.4) 

If the potential V(p) has degree larger than three, we obtain after a local 
expansion and a rescaling of fields, 

P-Ps = ( Qr ,l )(A - m 2 ) ex <p/N, x ^ Nx , (6.5) 

a simple super-renormalizable (p 4 (x) field theory. If we insist instead that 
the initial theory should be renormalizable, then we remain with only one 
candidate, the renormalizable (</> 2 ) 3 field theory, also relevant for the tricritical 
phase transition with O(N) symmetry breaking. Inspection of the potential 
W(p) immediately shows a remarkable feature: Because the term added to 
V(p) is itself a polynomial of degree three, the critical conditions lead to a 
potential W((p) which vanishes identically. This result reflects the property 
that the two saddle point equations (dS/dp = , OS/dX = in Eqs. (4.3)) are 
proportional and thus have a continuous one-parameter family of solutions. 
This results in a flat effective potential for (p(x). The effective action for tp 
depends only on the derivatives of ip, like in the 0(2) non-linear a model. 
We conclude that no non-trivial double scaling limit can be obtained* under 
these conditions. In three dimensions with a (</> 2 ) 3 interaction we can generate 
at most a normal critical point n = 3, but then a simple (<j) 2 ) 2 field theory 
suffices. 

The ambiguity of the sign of a(d) discussed in section 4 and in Appendix 
A2 has an interesting appearance in d = 3 in the small m 2 region. If one 
keeps the extra term proportional to a(d) in Eq. (6.3) we have 

w(p) = v(p) + 8 -^(pc- P f + a ^\ Pc - p)\ 

Using now Eq. (6.2) and, as mentioned in section 5, the fact that in the small 
m 2 region the potential is proportional to (p — p c ) 3 we can solve for m 2 . Since 



* It was pointed out that another sort of obstacle is present for (4> 2 ) 2 model 
in d = 2 [13] and d = A [14]. 
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m 2 > the appearance of a phase with small mass depends on the sign of 
a(d). Clearly this shows a non-commutativity of the limits of m 2 /A 2 — > and 
N — > oo. The small m 2 phase can be reached by a special tuning and cannot 
be reached with an improper sign of a(d). Calculated in this way, m 2 can be 
made proportional to the deviation of the coefficient of p s in V(p) from its 
critical value 16n 2 . 

7 Conclusions 

This is a study of several subtleties in the phase structure of O(N) vector 
models around multicritical points of odd and even orders. One of the main 
topics is the understanding of the multicritical behavior of these models at 
their critical dimensions and the effective field theory of the O(N) -singlet 
bound state obtained in the iV — > oo, g — > g c correlated limit. It is pointed 
out (in contrast to previous studies) that the integration over massive O(N) 
singlet modes is essential in order to extract the correct effective field theory 
of the small mass scalar excitation. After performing this integration, it has 
been established here that the double scaling limit of (</> 2 ) n vector model in 
its critical dimension d = 2n/(n — l) can result in a theory of a free massless 
O(N) singlet bound state. This fact is a consequence of the existence of 
flat directions at the scale invariant multicritical point in the effective action. 
In contrast to the case d < 2n/(n — 1) where IR singularities compensate 
powers of 1 /N in the double scaling limit, at d = 2n/(n — 1) there is no such 
compensation and only a noninteracting effective field theory of the massless 
bound state is left. 

Another interesting issue in this study is the ambiguity of the sign of a(d). 
The coefficient of m 2 A d_4 denoted by a(d) in the expansion of the gap equa- 
tion in Eqs. (4.3c) and (4.9) seems to have a surprisingly important role in 
the approach to the continuum limit (A 2 3> m 2 ). The existence of an IR 
fixed point at g ~ 0(A^ _1 ), as seen in the (3 function for the unrenormalized 
coupling constant in Eq. (4.15), depends on the sign of a(d). Moreover, as 
seen in section 3.1 the existence of a phase with a small mass m for the O(N) 
vector quanta and a massless O(N) scalar depends also on the sign of a(d). It 
may very well be that the importance of the sign of a(d) is a mere reflection 
of the limited coupling constant space used to described the model. This is 
left here as an open question that deserves a detailed renormalization group 
or lattice simulation study in the future. 
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Appendices 

Al Double scaling limit: d = 

We consider the sum Z of connected Feynman diagrams in d = dimension: 

n N / 2 r dx 



e z 



[ d N $e~ NV ^) = f te e _ N[v{x) _ ln{x)/2] ^ 

J T(N/2) J x 



where x is normalized by V(x) = \x + 0(x 2 ). 
We define 

W(x) = V(x) - iln(x). 
The saddle point equation reads 

W'(x s ) = 0. 

A critical point of order n is defined by the conditions: 

W{x)-W{x s ) ~ W (n \x s ) ( X - f s) 

The critical potential V n {x) of lowest degree thus satisfies 
V^x)=[l-(l-x/(n-l)r- 1 ]/(2x), 



and therefore 



px/(n-l) j 

Jo y 



and 



K ~ -2T^tt(i-^-i)) 



n-1 



x^n 2^=11 
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If we now change the variable, setting 



(l-x/in-l^N 1 / 71 = z, (A1.2) 

we find 



e Zc ~ 



r N/2 N -l/n 



T(N/2) exp[IJV(ln(n-l)-l^^ 

Adding to the critical potential all relevant perturbations, and shifting the 
saddle point back to x = n — 1, we obtain 



n-1 

V(x) = V n (x) + H [(1 " x/(n - l)) k - |(1 - s/(n - 1))»] . (Al.3) 
fc=i 

After the change of variable (Al.2), we find at leading order for N large 

n-1 

N{W(x) - W(n)) = ^ n + ^i - • 

fc=i 

We see that a double scaling limit is reached only when the coefficients Vk 
go to zero for large N with i> fc A 1_fc / n = Uk fixed. This implies that the 
coupling constants, i.e. the coefficients of V{x) in the expansion in powers of 
x must approach the critical values corresponding to the potential V n with a 
well-defined behaviour as a function of N. Note also that since for all Vk = 
several solutions of the saddle point equation coalesce, the critical potential 
corresponds to a singularity of the large N partition function in the space of 
coupling constants. Actually if we expand Z in powers of 1/N, 

Z = Y J N 1 ~ h Z h , 

h=i 

all terms Zh are singular at this point. This can be most easily seen by 
keeping only the most relevant term proportional to v\. Then one finds 

The scaling behaviour of v± as a function of N uses this singular behaviour 
to compensate the powers of N. Finally by tuning the coupling constants to 
reach a singularity we approach the radius of convergence of the perturbative 
expansion at h fixed, and thus enhance high order Feynman diagrams. Thus 
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one reaches a dense configuration of contributing Feynman graphs which is 
the analog of the sum over surfaces in O(N) matrix models. In the O(N) 
vector model this is an expansion over 'randomly branched polymers'. The 
scaling partition function is obtained by simplifying the integrand 

n-2 

N(W(x) - W(n - 1)) ~ + 1^ 

k=i 

because the Vk goes to zero and the term k = n — 1 can be eliminated by 
shifting z. Finally, if we keep only the most relevant perturbation term pro- 
portional to vi the double scaled partition function is given by a generalized 
Airy function 

Z( Ul )= I dzexp (-!(£+«!*)). 



A2 Regular izat ion and sign of the /3-function 

Here we give a few explicit examples which show the regularization depen- 
dence of the non-universal constant a(d) defined in (4.9). 

Pauli-Villars regularization. For the regularization (3.2) it is given by 




Let us consider for example 

D(k 2 ) = l + ak 2 + [3k 4 . 

D(k 2 ) > implies 4(3 > a 2 . 

If a > then clearly for every k, D(k 2 ) > 1, and a(d) is positive. However 
if a < it is possible to choose the parameters in such a way that a(d) will 
change sign; especially if one takes a close to — 2-/]3. Possible values of a(d) 
are exhibited in Fig. 1 below. 
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d 



Fig. 1: Values of a(d) in Pauli-Villars regularization. 



(From top to bottom: a = 0.8, 0.4, 0, -0.4, -0.8, -1.2, -1.6 and (3 = 1.) 

Lattice regularization. We extend here the calculations of [12] to more 
general lattice regularizations, to understand how general this sign property 
is. 

An alternative representation to Eq. (4.4) is: 



where (Q(k)) 1 is the massless free field propagator. Using the asymptotics 




of f(a) 



/(«) 



a 



oo 



(27m)-2 [1 + 0(1)], 



then 



I 



2/2 -l) (/(a) -(2™)"^). 



P ~ Pc = - 



K(d)m d - 2 + \ 



da ( e 



— am 



Therefore a(d) in Eq.(4.9) reads: 



a 



poo 

(d) = -\ daa(f(a)-(27ia)- d/2 
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a(d) depends on the regularization used for the lattice action 



s($) = s k ($) + nJ2[v($ 2 )] 



where 



S K ($) = \N{2n) d I d d k$(-k)Q(k)$(k), ${x) = j d d k e lkx $(k). 

J— 7T J— 77 



The lattice spacing is taken to be 1. 

If only nearest-neighbor interactions are included: 



then 



Q(fc)=4^sin 2 (%), 



f(a) = (e-«I (a)y 



where Io(ot) is a modified Bessel function. 

In the following regularization which includes contributions from next to 
nearest neighbors and next to next to nearest neighbors the dimension d can 
be varied continuously: 



(l-r)£ (0(z + A) + <t>{x -p.)- 2</>(x)) 



a.=i 



+ ±r (<t>{x + 2 ft) + <p(x - 2ft) - 2cf>(x)) 



The inverse of the free field propagator is now: 



a 

Q(k; r) = ( 4 sin2 ¥ - 4r sin4 • 



ai=i 



Values of a(<f) in this regularization are exhibited in Fig. 2 below. 
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Fig. 2: Values of a(d) in lattice regularization. 
(From top to bottom: r = -1.2, -0.8, -0.4, 0, 0.4.) 
At d = 3, including all next-to-nearest-neighbors we have: 



(i-OE + A) + 00* - A) - 20(x)) 



-jr ^ (0(x+/i+i>) + 0(x+/i — + /t+i>) + 0(x— /i— i>)— 40(x)) 

3>A*>^>1 



and thus 



Q 2 (£;;r)=4^sin 2 (%)-2r 



Once again, one finds that as the parameter r is growing, a(d) is changing its 
sign from a(d) > to a(d) < ( we have a(d) = 0.002 , -0.0001 , -0.07 for 
r = —0.6 , —0.5 , 0.9 respectively ). 
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